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oe binary variable indicates on/off situation of DG units 
er al binary variable indicates charge/discharge situation of 
ES units 
SOE wi state of the energy level of ESs (kWh) 


5.1 INTRODUCTION 


In the context of future smart grids, the distributed and small-scale energy resources 
will play a more remarkable role in comparison with their present situation. One of 
the main technical and commercial challenges for integration of distributed energy 
resources (DERs) into the grid is their invisibility from system operator’s viewpoint. 
Virtual power plant (VPP) concept seems an appropriate solution to tackle such a 
problem. In fact, a VPP is composed of small DERs (traditional and renewables), 
controllable loads, and energy storage (ES) devices and can coordinate them in order 
to provide several local and system services through participation in different elec- 
tricity markets with the aim of maximizing its profits. 


5.2 OBJECTIVE 


As mentioned before, the VPP aggregates and manages not only distributed genera- 
tions such as diesel generators (DGs), photovoltaic (PVs), and wind turbines (WTs) 
but also ESs as well as flexible loads (FLs) through an integrated optimization 
framework so that there is a maximization of its revenues in various electricity mar- 
kets. In other words, a VPP is a flexible representation of a portfolio of DERs that can 
be used to make contracts in the wholesale market and offer services to the system 
operator. On this basis, the objective of this deliverable is to develop a centralized 
optimization engine for a VPP that can be operated in different modes, including 
technical VPP (TVPP) mode, commercial VPP (CVPP) mode, and also energy com- 
munity (EC) mode according to the user preferences. 


5.3. MATHEMATICAL MODEL OF OPTIMIZATION 
ENGINE IN DIFFERENT OPERATION MODES 


5.3.1 TVPP MopeL 


A TVPP is a type of VPP that consists of DERs from the same geographic location. 
The TVPP encompasses the real-time impacts of the local network on the aggre- 
gated profile of DERs as well as representing the cost and operating constraints 
of the portfolio. The TVPP can participate in day-ahead and intraday energy mar- 
kets for purchasing or selling power on behalf of the existing resources under its 
umbrella. The operator of a TVPP, so-called VPP aggregator, requires detailed 
information of the local network as well as the existing resources [1]. In fact, in 
this case, the optimization engine must be considered the network structure and 
run an optimal power flow in order to avoid either feeder’s congestion or voltage 
magnitude violation. 
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To model the stochastic nature of renewable resources as well as taking into 
account the electricity market price’s volatility, the proposed methodology makes 
use of a two-stage stochastic programming framework and considers a network- 
constrained market participation procedure from a TVPP point of view. The struc- 
turing of the stochastic programming problem into two stages is justified by the fact 
that the TVPP aggregator must submit its offers into different electricity markets, 
while facing a number of uncertainties such as volatility in market prices, load con- 
sumption, and generation of renewable resources. The proposed model has a hierar- 
chical structure so that it encompasses both day-ahead and intraday time frames as 
can be observed in Figure 5.1. 

In fact, the optimal strategy of VPP for participation in day-ahead energy market 
is determined at the first stage, while the second stage is assigned for intraday market 
decisions which are based on the updated information obtained between the closure 
of the day-ahead and the intraday energy markets. 

The objective function of the TVPP is formulated in Eq. (5.1). The objective func- 
tion is the total expected cost during the scheduling horizon, which contains three 
terms, including cost/revenue in day-ahead energy market, cost of DG power pro- 
duction, and incentive payment to FL, respectively. It is notable that the objective 
function that must be optimized for the second stage is completely similar to Eq. (5.1) 
with the difference that the power traded in day-ahead market is an input parameter 
since the day-ahead market has been cleared. In fact, the first term in the objective 
function must be substituted by the revenue or cost term for selling or purchasing 
power in intraday market. 
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FIGURE 5.1 Schematic of the proposed hierarchical optimization model. 
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The active and reactive power balance constraints have been formulated in Eqs. (5.2) 
and (5.3), respectively. Note that the first term in Eq. (5.2) is associated with power trans- 
actions with the upstream grid and just must be taken into account for the source bus. 
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The active and reactive power relations are considered according to the power factor 
concept as can be observed in the following equation. 


Ore a PF x PP vt 
Ost, = PF x Phi (5.4) 
OF = PF RPE 


The TVPP must consider the power flow restrictions within its optimization model. 
On this basis, the linear form of load flow equations has been formulated in Eqs. 
(5.5)—(5.15) [2]. 
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Equation (5.5) is considered with the aim of balancing voltage between two nodes. 
It is notable that V2 in Eq. (5.5) is an auxiliary variable that shows the linear form 
of squared voltage relation. Also, linearization of active and reactive power flows 
that appear in the apparent power is formulated in Eq. (5.6). Equations (5.7)—(5.11) 
have been formulated for piecewise linearization. The number of blocks needed to 
linearize the quadratic curve is considered to be five based on [3], which maintains 
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the right balance between accuracy and computational requirements. Further 
descriptions and justifications of the network model used can be found in [4, 5]. The 
maximum allowable current flow of branches is taken into account in Eq. (5.12). 
Note that /2 refers to an auxiliary variable that demonstrates linear form of the 
squared current flow F in a given branch. Moreover, at most one of these two positive 
auxiliary variables, i.e., P, w: and Q; p w1 can be nonzero in a time. This condition 
is again implicitly enforced by optimality. Equations (5.13) and (5.14) restrict these 
variables through the maximum apparent power for completeness. Finally, Eq. (5.15) 
represents the allowable voltage magnitude at each node. The remaining constraints 
are associated with the existing resources within the TVPP, which are presented in 
the following. 


5.3.2 CVPP MopeL 


CVPP performs commercial aggregation and does not take into account any network 
operation aspects that active distribution networks have to consider for stable opera- 
tion [6]. The aggregated DER units are not necessarily constrained by location but 
can be distributed throughout different distribution grids. Hence, a single distribu- 
tion network region may have more than one CVPP-aggregating DER units in its 
region. The objective function of CVPP is completely similar to Eq. (5.1); however, 
the active power balance constraint is different due to non-consideration of network. 
Also, there is no need to consider reactive power in the case of CVPP due to the same 
reason. The active power balance constraint for CVPP is modelled in Eq. (5.16). 


DA DG PV Wind Dch Ch les D 
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į pv wf es fl 
The remaining constraints are assigned to the existing resources within the CVPP, 
which are presented in the following. 


5.3.3 ENERGY COMMUNITY MODEL 


EC is an innovative concept that has been developed to facilitate grid integration 
of small-scale renewable energy resources, optimize the consumption pattern of 
customers, and alleviate the loading of the grid through using available flexibility 
of active prosumers. In the following, a general model of an EC consisting of sev- 
eral small-scale resources is given. It is notable that these resources are spatially 
very close to each other and connected to the same distribution network. Also, it is 
assumed that there is a non-profit community coordinator who dispatches the exist- 
ing resources within the community and manages the transactions with upstream 
utility as well as other ECs in order to supply the load with minimum cost. The con- 
ceptual schematic of the considered model is presented in Figure 5.2. 

As it can be observed in Figure 5.2, the EC coordinator runs the local market in 
order to find an optimal dispatch for the components within the EC, including PV 
generation, FLs, DGs, and ES units. In this regard, the ECs can also have transac- 
tions with other ECs in order to share their available flexibility with each other. 
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FIGURE 5.2 Conceptual schematic of energy community model. 


Moreover, the ECs can buy or sell power to the upstream utility with the aim of 
ensuring power balance within the community. The mathematical formulation of 
EC is presented later. 

The objective function is a cost minimization problem that has been formulated 
in Eq. (5.17). The first and second terms in the first line of the objective function 
represent the costs/revenues as a result of transactions with the upstream utility. The 
first term in the second line of Eq. (5.17) associates with the cost of DG’s power gen- 
eration containing their start-up cost. Finally, the last term of objective function is 
assigned to incentive payments for load curtailment through FLs. 
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The A in the previous equation represents the time step coefficient and can be defined 
as A = Time _ step/60. The Time_step in fact indicates the resolution of optimization 
time slots in minutes. The previous objective function should be minimized while 
satisfying a number of constraints for different devices. The EC power balance con- 
straint has been modelled in Eq. (5.18). 
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5.4 MATHEMATICAL MODEL OF COMPONENTS WITHIN THE VPP 


This section is assigned to the model of existing resources within the VPP and is 
similar for all operation modes, i.e., TVPP, CVPP, or EC. The considered compo- 
nents are DGs, ESs, and FLs that are modelled in this section. Note that the renew- 
able generations such as PV and WT have been taken into account. According to the 
priority of renewable resources in dispatch, these resources are modelled to such a 
negative demand. In fact, the whole production of these resources is integrated into 
the VPP power scheduling. Moreover, the operation costs of the mentioned resources 
assume to be zero that is a logical assumption. 


5.4.1 DG Units 


The constraints in Eqs. (5.19)—(S.25) are related to technical restrictions of DGs. The 
minimum and maximum ranges of output power of DGs are shown in Eq. (5.19). 
The ramp-up and ramp-down capabilities of conventional DGs are formulated in 
Eqs. (5.20) and (5.21), separately. Equation (5.22) is assigned to start-up cost of DGs. 
Moreover, Eqs. (5.23) and (5.24) indicate minimum up and down time limitations of 
DGs, respectively. 
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5.4.2 ES Units 


The constraints of ESs that should be satisfied are formulated in Eqs. (5.25)—(5.30). 
Limitations on charging and discharging power of ESs are modelled in Eqs. (5.25) 
and (5.26). The state of energy (SOE) level and its limits are also given in Eqs. (5.27) 
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and (5.28), separately. Furthermore, it is presumed that the SOE level at the end of 
scheduling period must be greater or equal to the amount of SOE at the beginning of 
scheduling period as formulated in Eq. (5.29). This is due to the fact that the ES units 
must have sufficient energy for the next day. 
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5.4.3 FLEXIBLE LOADS 


The other set of constraints is associated with FLs as shown in Eqs. (5.30)—(G.32). 
Equation (5.30) restricts the amount of FL due to the fact that just a portion of the 
load is flexible. Also, the maximum load-pickup and load-drop rates for flexible 
demand are stated in Eqs. (5.31) and (5.32), respectively. 
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5.5 NUMERICAL ANALYSIS AND DISCUSSIONS 


In order to show the effectiveness of the proposed model, some numerical analyses 
have been reported in this section. To this end, a CVPP is considered with the fol- 
lowing assets. There is a 500-kW wind farm and a 200-kW PV site so that each of 
the wind farm and PV site generation is modelled considering three independent 
scenarios, including as forecast, high, and low, with probabilities 0.6, 0.2, and 0.2, 
respectively. Also, there is one ES unit with the energy capacity of 200 kWh and 
maximum charging/discharging rates of 100 kW/h. The charge and discharge effi- 
ciencies of ES are assumed to be 85%. In addition, it is assumed that the maximum 
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TABLE 5.1 
DGs Technical and Cost Data 


DG 1 600 100 105 120 20 0.040 
DG2 950 200 175 200 30 0.035 


and minimum SOE of the ES are equal to 90% and 10% of its energy capacity due to 
life-time considerations. The initial SOE of ES is 50% of its energy capacity. 

It is assumed that only there are three FLs so that they can curtail 20% of their ini- 
tial consumption in each hour in response to 0.035 $/kWh as an incentive payment. 
The load-pickup and load-drop rates are considered to be 25 kW/h. Furthermore, two 
conventional DGs are available with technical characteristics as well as cost terms 
reported in Table 5.1. 

Figure 5.3 represents the trading power of CVPP in the day-ahead market. In 
this figure, the positive values are assigned to purchased power, while the negative 
values represent the sold power in the day-ahead market. According to this figure, 
it is obvious that the CVPP purchases power in time periods with low market prices 
(early morning and end of the night), whereas it sells power to the market during the 
day when the market price is relatively high. 

The obtained results for the optimal dispatch of various components belonging 
to CVPP have been reported in Figure 5.4. According to Figure 5.4, when the aggre- 
gated generation is more than consumption, the CVPP sells the extra power to the 
market and vice versa. It is obvious that both DGs are generating so that their genera- 
tion in the early morning and end of the day is set to minimum, and during the day, 
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FIGURE 5.3 CVPP trading power vs. expected day-ahead market prices. 
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FIGURE 5.4 Power dispatch of components belongs to CVPP. 


these DGs generate with full capacity due to the fact that their generation costs are 
lower than market prices. Moreover, it can be seen that the CVPP utilizes the FLs at 
all the periods. 

The charging and discharging plan of ES during the scheduling horizon has been 
illustrated in Figure 5.5. As it can be seen, the ES unit charges during the early 
morning period, particularly at hours 4:00 and 5:00, and injects its power to the grid 
during the day, specifically between 7:00 and 11:00. The charging and discharging 
plan of ES unit is completely based on the day-ahead market variations. In fact, in 
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FIGURE 5.5 Charging power, discharging power, and state of energy of ES. 
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FIGURE 5.6 Power dispatch of components belongs to CVPP with emission cost 
consideration. 


addition to ES role in facilitating the integration of renewable generations such as PV 
and WT, it uses the energy price arbitrage and makes more profit for CVPP. 

The other important points that can affect the CVPP strategy is consideration of 
pollutant emission cost for DGs due to the fact that the DG’s technology is mainly 
based on diesel burning. In this case, it is assumed that the costs of DGs are multi- 
plied by 1.3. On this basis, both of the DGs never generate and the CVPP is forced 
to purchase power from market at all time. The power dispatch in this case is shown 
in Figure 5.6. 

As mentioned in the modelling part, the stochastic programming approach has 
been used here in order to take into account the uncertainties of renewables, load, 
and market price. In order to evaluate the impacts of forecasted parameters on the 
power dispatch of various resources, the obtained results for the real case are com- 
pared with outputs of the stochastic model. The forecasted scenarios for PV and WT 
have been compared with the real generation in Figure 5.7. 

In such a situation, the obtained results for the cost terms of objective function 
have been compared in Table 5.2. 


TABLE 5.2 
Comparison of Cost Terms in Objective Function 
Market Transaction Flexible Load 
Case Cost ($) DG Cost ($) Cost ($) Total Cost ($) 
Scenario-based —586.5 965.6 420.7 799.8 


Real data —588.9 973.1 414.0 798.2 
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FIGURE 5.7 Scenarios and real data of renewable generation: (a) Wind generation (p.u.) 
and (b) PV generation (p.u.) 


The obtained results for market transactions of the CVPP in two mentioned cases 
(real and forecasted data) have been compared in Figure 5.8. As observed, the pro- 
posed model is relatively robust in the face of forecast tool’s error. However, there are 
some deviations in the number of time slots (red box). 
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